New Measurements and Quantitative Analysis of Electron Backscattering in the 

Energy Range of Neutron /9-Decay 
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We report on the first detailed measurements of electron backscattering from plastic scintillator 
targets, extending our previous work on beryllium and silicon targets. The scintillator experiment 
posed several additional experimental challenges associated with charging of the scintillator target, 
and those challenges are addressed in detail. In addition, we quantitatively compare the energy 
and angular distributions of this data, and our previous data, with electron transport simulations 
based on the Geant4 and Penelope Monte Carlo simulation codes. The Penelope simulation is found 
globally to give a superior description of the data. Such information is crucial for a broad array of 
weak-interaction physics experiments, where electron backscattering can give rise to the dominant 
detector-related systematic uncertainty. 

PACS numbers: 34.80.Bm, 23.40.-s, 29.30.Dn 



I. INTRODUCTION 

Our previous work on backscattering from beryl- 
lium and silicon extended the work of others 0, 0, El to 
energies relevant for nuclear physics applications. We fur- 
ther extend this work to a more relevant target, organic 
scintillator. The main experimental challenge arose from 
the fact that the target material is not conducting, ne- 
cessitating mitigation of effects due to charging in order 
to accurately sense currents. 



II. EXPERIMENT OVERVIEW 

As for our previous measurements, the experiment 
consisted of an electron gun and a scattering chamber 
containing a movable target Two modes of acquir- 
ing backscattering data were used. In one mode a sil- 
icon detector was used to detect the energy and angle 
of backscattered electrons. In a second higher-current 
mode, the electrical current due to the backscattered elec- 
trons incident on the chamber walls was measured. These 
two modes were referred to as silicon detector mode and 
current integration mode, respectively. 

The plastic scintillator targets used for the experiment 
were obtained from Eljen Technologies 0- The type of 
plastic scintillator used was EJ-204. The density of EJ- 
204 is 1.032 g/cm 3 , with the dominant elemental compo- 
sition being 5.21 x 10 22 H atoms/cm 3 , and 4.74 x 10 22 C 
atoms/cm 3 . A sample of plastic scintillator was coated 
with a 500 A thick layer of aluminum via evaporation. 
Targets of a suitable size and shape were then cut from 
this sample, with care taken to preserve the aluminizcd 
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front face. The resulting targets were l"xl" square, 
3 mm thick, with a thinner tab of aluminized plastic 
projecting from the top for suspension from a central 
rotatable feedthrough. 

As before, the energy and angular distribution of the 
backscattered electrons was measured using a 25 mm 2 
active area silicon detector. The energy resolution of the 
silicon detector system was typically a — 4.3 keV, and 
was independent of energy. 



III. BACKSCATTERING MEASUREMENTS 

Backscattering measurements were performed for nor- 
mal incidence upon the target. Measurements were per- 
formed for incident electron energies of 43.5, 63.9, 83.8, 
104, and 124 keV. For each energy, both silicon detector 
mode and current integration mode measurements were 
taken. 



A. Silicon Detector Mode 

The quantities which varied for the silicon detector 
measurements conducted at a particular electron beam 
energy .Ebcam were the detection angle 8, and the energy 
E the backscattered electron deposited in the silicon de- 
tector. The angle 9 was defined relative to the surface 
normal vector of the material, so that 8 — degrees corre- 
sponds to a backscattering event where the electron went 
directly back along the incident beam direction. The di- 
mensionless energy q = E/Eb eaxn will also be used. 

For silicon detector mode, the systematic uncertainty 
is unchanged from our previous work, resulting in an av- 
erage 12% normalization systematic for each angular set- 
ting of the detector. However, we specifically addressed 
current detection and charging of the target, as will now 
be discussed. 
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To ensure the reliability of measurements of current, 
we compared the measurements of two different, care- 
fully calibrated current integrators against one another 
for identical experimental conditions, and in turn cross- 
compared those measurements against two calibrated pi- 
coammeters. The measurements were consistent at the 
3% level. 

The possible effects of charging and incomplete current 
detection were monitored by observing scintillation light 
from the electron beam as it struck the scintillator target, 
on a camera mounted behind the target outside a view- 
port in the chamber (at 9 — 180°). Over the course of 
taking a complete angular range of data for a particular 
incident beam energy, the brightness of the scintillation 
light did not visibly change with time. Previous tests 
with uncoated scintillator, or scintillator coated poorly 
with graphite, had shown that the emanation of scintilla- 
tion light from the target would eventually cease, indicat- 
ing that the electron beam had been diverted away from 
the central spot on the target. Also discharges would 
be seen due to arcing from the face of the scintillator to 
the conducting target rod. No such effects could be seen 
with the Al-coated scintillator target. On leaving the 
beam on the target for long periods of time, no change 
in the sensed current was seen at the 1% level. To also 
search for charging, the electron beam could be switched 
off and on rapidly by inserting a Faraday cup upstream 
of the chamber. Upon restoration of the beam, the cur- 
rent sensed by the target was found to agree with the 
value before intercepting the beam with the Faraday cup 
to the 1% level. 

The effect of the aluminum layer can be estimated by 
comparing the thickness of the Al layer to the mean range 
of the electrons. The mean range for 43.5 keV electrons 
in plastic scintillator is 30 /im p , which is three orders of 
magnitude larger than the thickness of the Al layer. This 
implies that scattering from the Al is suppressed below 
the 1% level, even for the lowest energy incident electrons 
reported in this work. Both estimates were confirmed in 
Monte Carlo studies of scattering from thin layers 0. 
For the silicon detector data, the effects of the Al layer 
were therefore neglected. 

As a cross-check of the reproducibility, dead-time, and 
current detection uncertainties due to charging, runs 
were taken for various beam currents. The normalized 
yield for these runs was found not to vary outside the 
previously quoted 7% reproducibility uncertainty. 

For completeness, a catalog of the dominant system- 
atic uncertainties considered for silicon detector mode is 
displayed in Table [I] Aside from target deterioration, 
which will be discussed in the next section, these sys- 
tematic uncertainties were described in detail in our pre- 
vious work 0| . The effects listed in the upper portion of 
Table prefer to results for the observable -^-J^, and 
average 12%. The effects in the lower portion of Table |I| 
refer to extrapolation uncertainties encountered when in- 
tegrating these data over q and/or 6, and these must be 
added in quadrature to those above when considering our 



Effect 


Uncertainty 


Reproducibility 


7% 


Active Area 


4% 


Finite Beam Spot 


5%x sin6> 


Dead Time 


3% 


Alignment 


2% 


Current Detection 


3% 


Target Deterioration 


1% 


Extrapolation over q 


6%-20% 


Extrapolation over 9 


4% 


Total 


12-23% 



TABLE I: Dominant systematic uncertainties for silicon de- 
tector mode. 

results for drj/d£L and ry from silicon detector mode. For 
extrapolation over q, the uncertainty is larger for lower 
beam energy, as the finite detection threshold becomes 
more important; hence the 20% value refers to beryllium 
at 43.5 keV incident energy, since the beryllium data is 
more peaked at lower q, while the 6% value refers to 
data taken with 124 keV incident energy. Considering 
the range given for the total systematic uncertainty, the 
lower bound (12%) refers to all observations of , 
while the upper bound (23%) refers specifically to the 
extraction of r\ from beryllium at 43.5 keV, where ex- 
trapolation uncertainties dominate. 

For these experiments, the silicon detector threshold 
varied from 10 keV to 18 keV depending on noise condi- 
tions at the time of any given data-taking run. 

Effects due to backgrounds from X-rays, and multiply 
backscattered electrons were also discussed in detail in 
our previous work. Multiply backscattered electrons re- 
sulted in systematic uncertainties at the 3-5% level at 
q = 0.2 for the highest beam energy considered. But this 
contribution to the systematic uncertainty is below 1% 
above q — 0.3, and is negligible for lower beam energies. 



B. Current Integration Mode 

For each beam energy setting, current integration 
mode measurements were also performed. 

The largest potential systematic effect, which was pre- 
viously uncontrolled, arose from deterioration of the scin- 
tillator target by the electron beam. Over the course of 
an hour, running at beam currents of tens of nA, r\ was 
observed to steadily increase with time, plateauing at a 
value typically 15% larger than its initial value. The tran- 
sition to larger r\ was found to occur more rapidly with 
higher beam currents and higher beam energies. The 
glow of the beam spot on the target was also monitored 
on the video camera and found to decrease in bright- 
ness in a correlated way. The brightness was found not 
to recover after leaving the scintillator in vacuum over 
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several days, ruling out charging/discharging of the bulk 
scintillator. Upon removal of an affected target from the 
vacuum chamber, a small brownish spot within the scin- 
tillator could be observed, with no obvious deterioration 
of the mirror-like aluminized front face. We believe the 
discoloration, reduction in scintillation light, and increase 
in i] are symptoms of a chemical change in the scintilla- 
tor, possibly resulting in the liberation of hydrogen. This 
would increase the relative carbon content and hence r]. 

To control this potential systematic, new scintillator 
targets were used and were exposed to electron beams 
with currents less than 1 nA. The targets were exposed 
to beam for the minimum time possible for currents to be 
sensed accurately using picoammcters. This limited the 
contribution to the systematic uncertainty to less than 
1%, based on the slope measured during the prior hour- 
long measurements. 

A dependence on beam current at the 3% level had 
been seen for our previous current integration mode data 
on Be and Si target jl]. For those data, the current de- 
pendence was attributed to charging of various compo- 
nents in the chamber. As our new studies were done 
generally at lower beam currents, we expect this contri- 
bution to be smaller, but retain the pessimistic upper 
limit of 3%. 

In order to characterize and control the effects of low- 
energy secondary electrons (for our work, defined to be 
below 50 eV in energy), a wire mesh cage was inserted 
into the chamber and held at -55 V. From an electrostatic 
model of the electric potential in the chamber, this re- 
sulted in a minimum potential wall of -50 V between the 
target and the chamber. For the purposes of modelling 
these measurements, the integrated backscatter current 
was therefore taken to be due to all electrons emanating 
from the target with greater than 50 eV. 

The wire mesh cage was referred to as "the grid" con- 
sistent with terminology used in previous backscattering 
literature. The grid was constructed from a thin copper 
rod bent into a cylindrical shape. Steel wire was wound 
on the resultant frame in an end-over-end pattern result- 
ing in vertical wires equally spaced running down the 
sides of the cylindrical shape. The resultant cylindrical 
wire cage resembled a bird cage which enclosed the tar- 
get. In this way, biasing the grid at negative potential 
prevented secondaries from traveling from the chamber 
walls to the target and vice versa. 

However, as with our previous work, the insertion of 
the grid resulted in systematic effects due to incom- 
plete compensation for secondaries and due to effects of 
backscattered electrons striking the grid itself. 

A residual correction due to a piece of conducting tar- 
get rod penetrating the top of the grid must be applied. 
The rod subtended a small but finite solid angle viewing 
the beam spot on the target. High-energy backscattered 
electrons could strike that portion of the rod, resulting 
in secondary electron production. Secondary electrons 
created on that portion of the rod would not be sup- 
pressed by the grid and would be collected on the target. 



This resulted in a residual dependence of the apparent 
1] on target voltage, which could be corrected. The cor- 
rection gave rise to a 7% contribution to the systematic 
uncertainty in the previous work, and was the dominant 
uncertainty. For this work, the contribution was reduced 
to 3% by reducing the solid angle subtended by the rel- 
evant portion of the target rod. This reduced the size 
of the correction, and hence the systematic uncertainty. 
Detailed analytical calculations of the effect of the solid 
angle and its variation with e.g. distance of the beam 
spot from the top of the target close to the target rod 
gave additional confidence in this uncertainty. As be- 
fore, secondaries created on the grid were characterized 
by monitoring the grid current, and contributed at the 
1% level. 

A potential systematic effect arises due to differences 
in backscatter yields above 50 eV because of the presence 
of the Al coating on the target. Electrons with energies 
above 100 eV could not be adequately assessed via al- 
tering the voltage on the grid, while from the electron 
transport arguments presented earlier, only for energies 
above about 10 keV can the effects of the Al layer be 
argued to be negligible. From similar considerations to 
the extrapolation uncertainty for silicon detector mode, 
we limit possible effects of the Al coating in this energy 
range to the 6% level. 

The dominant systematic uncertainties for current- 
integration mode results for our new scintillator target 
data are listed in Table ITT1 



Effect 


Uncertainty 


Target rod correction 


3% 


Grid secondaries 


1% 


Reproducibility 


5% 


Current dependence 


3% 


Target deterioration 


1% 


Al coating effects 


6% 


Total 


9% 



TABLE II: Dominant systematic uncertainties for current in- 
tegration mode for scintillator target data. 



IV. RESULTS 
A. Silicon Detector Mode 

The normalized, background-subtracted spectra accu- 
mulated for various detector angles for 124 keV electrons 
normally incident on a scintillator target is shown in 
Fig. ^ The data are plotted as a function of the di- 
mensionless energy q. On the vertical axis, jj- f q ^ Q , the 
number of counts per incident electron, per unit q, per 
unit solid angle is plotted. In the absence of the effects 
of detector response (resolution and backscattering), this 
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FIG. 1: (Color online) Normal incidence backscattering from 
scintillator target at Ekeam = 124 keV. Curves represent data 
taken with silicon detector. Histogram is Monte Carlo sim- 
ulation based on (a) Geant4, and (b) Penelope. Systematic 
uncertainty in the normalization of the data is estimated to 
be 12% on average, ranging from 11% at small angles to 15% 
at large angles. A scale factor of 1.073 is applied to the (a) 
Geant4 simulation. In (b), the factor is 1.163 for Penelope. 



would be the normal-incidence backscattered fraction per 
unit q, per unit solid angle. 

Monte Carlo simulations of backscattering were car- 
ried out with custom codes based on the Geant4 y| and 
Penelope toolkits. The custom aspects of the simula- 
tion codes were described in Ref. for example, their 
handling of backscattering from the silicon detector it- 
self and the detector response function (relating to this 
particular observable). Fig. ^ compares the data with 
the results of these two codes. In each scale fac- 

tor is applied to the Monte Carlo. The scale factor was 
determined from a fit to the data, which is described in 
Section [Vj As for our previous work, Geant4 systemati- 
cally underestimates the peak in the data near q = 0.95. 
However the positions in q of the low-energy and elas- 
tic peaks are rather well-described by Geant4. In the 
case of the Penelope simulation, when the Monte Carlo 
is rescaled, it is apparent that trends in both energy and 
angle are well represented by Penelope. This result is 
discussed quantitatively in Section Ivl 

The angular dependence of the backscattered fraction 
can be determined by integrating the data over q. The 
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FIG. 2: (Color online) dr//dQ for beryllium (triangles) and sil- 
icon (inverted triangles) targets at i?beam = 124 keV. Black 
points with error bars indicate data with total normalization 
systematic uncertainties shown. Red solid histogram indi- 
cates the results of the Geant4-based Monte Carlo simula- 
tion. Green dot-dashed histogram indicates the results of 
the Penelope-based Monte Carlo simulation. No Monte Carlo 
scale factors are included. 



result of doing so is shown in Fig. [21 and is compared 
with our previous results for silicon and beryllium tar- 
gets. A linear fit based on the first 20 keV of data above 
the analysis cut was used to extrapolate to 50 eV (the de- 
fined threshold for secondary electrons), so that these in- 
tegrals and subsequent integrals could be compared with 
the current integration measurements. An additional sys- 
tematic uncertainty was assigned to the extrapolation, 
based on differences of fit functions, and comparison to 
models. For 124 keV beam energy, this extra systematic 
uncertainty was 5%, resulting in a contribution due to 
extrapolation of typically 6% (see Table E}. 

Fig.|2]also compares the Geant4 and Penelope simula- 
tions with the data. No Monte Carlo scale factors are ap- 
plied for this comparison. Each Monte Carlo separately 
correctly predicts that the scintillator results should be 
larger than the beryllium results, as expected due to the 
larger Z of the carbon nuclei in scintillator. As noted pre- 
viously, the Penelope simulation tends to better describe 
trends in angle. The Geant4 distributions are somewhat 
narrower compared to the data and Penelope. Addition- 
ally, the Geant4 simulation gives systematically larger 
backscattering from each material than does the Pene- 
lope simulation. 

The data were integrated over angle to determine the 
total normal-incidence backscattered fraction rj. The re- 
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suits under varying conditions (see Table ITT|) . In the case 
of the silicon detector measurements, it is dominated by 
reproducibility, and by uncertainties in alignment, solid 
angle effects, and extrapolation to 50 eV (see Tabic QJ. 

The data are also compared with previous data on Be 
and Si targets due to Drescher et al. [lij and with data on 
Be due to Neubert et al. [H|. Both groups used current 
integration techniques to arrive at their results. Neubert 
et al. |Tl| in particular used a second target apparatus 
to study the effects of secondary electrons, as opposed 
to the grid used in this work. Only the subsets of their 
data that overlap the region 43.5 to 124 keV are plotted. 
The Drescher data on Be are systematically higher than 
the Neubert data. As noted previously, our data tend 
to agree with the Neubert data on Be, as do the data of 
Massoumi et al. 0, 3 taken below 40 keV. 

Our data on organic plastic scintillator to our knowl- 
edge are the first in this energy range. As expected, the 
results lie below the previous measurements on carbon 
(not shown), and above our data and the Neubert data 
on beryllium. 



FIG. 3: (Color online) Normal incidence backscattering from 
Be, Si, and plastic scintillator targets. Integrated silicon de- 
tector measurements are shown by the inverted filled trian- 
gles. Current integration measurements are shown by filled 
circles. Total systematic uncertainties are shown and the cur- 
rent integration measurements are displaced by 2 keV so that 
the error bars do not overlap. Previous current integration 
measurements [Tol ITl| are displayed. The histograms show 
the results of the Geant4 (Red solid) and Penelope (Green 
dot-dashed) Monte Carlo simulations. 

suits of this integration are shown in Fig. [3| and are com- 
pared with current integration measurements (described 
in section ITV B|l . with previous data, and with the mod- 
els. In Fig. |3| the total systematic uncertainty, including 
extrapolation to 50 eV and extrapolation over unmea- 
sured angles, is plotted. 

Fig-Elcompares the data with the Geant4 and Penelope 
simulations. The same discrepancies in normalization are 
again observed. Both Penelope and Geant4 adequately 
describe the reduction of r\ as the beam energy increases. 

B. Current Integration Mode 

The results for r\ based on our current integration mea- 
surements are also shown in Fig. 

The silicon detector measurements are found to be sys- 
tematically higher than the current mode measurements; 
however, the two methods of are found to agree within 
the systematic uncertainties. In the case of the current 
integration method, this systematic uncertainty is dom- 
inated by residual correction for secondary electron col- 
lection due to the penetration through the grid of the 
target rod and reproducibility of the experimental re- 



V. QUANTITATIVE COMPARISON TO 
MODELS 

Monte Carlo to data fits including a single fit param- 
eter, a global normalization factor, were performed for 
four different observables measured by our experiments 
(including those presented in Ref. Q). The observables 
considered for fitting were: ^ d at 124 keV beam en- 
ergy, ^ at 124 keV beam energy, r\ from silicon detector 
mode, and rj from current integration mode. The tech- 
nique of x 2 minimization was used to constrain the fit. As 
mentioned in our previous work , point-to-point statis- 
tical and systematic uncertainties were exceedingly small 
relative to the overall normalization systematic. To sim- 
plify the fitting procedure, it was assumed that the point- 
to-point uncertainties were proportional to the global 
normalization uncertainty, for the purposes of evaluat- 
ing x 2 - 

The results of the fit are listed in Table E2 The fit 
results for d ^ (124 keV) data for scintillator are dis- 
played in Fig. Note that in our previous work 
we did not use this fitting technique, and normalization 
factors were simply determined by eye. However, the nor- 
malization factors determined using the x 2 minimization 
method have a good correspondence with those numbers. 

As the absolute point-to-point uncertainty was not de- 
termined precisely, the values of uncertainties on the 
fit parameter, and of absolute x 2 ' s ; have no mean- 
ing. Therefore only the ratio of % 2 's determined for the 
Geant4 model, divided by that for the Penelope model is 
quoted. I.e., when taking the point-to-point uncertainty 
to be equivalent to the total normalization uncertainty, 
the value of the reduced x 2 was generally significantly 
less than unity. The exception to this was the compar- 
ison of the Geant4 simulation to -A- J^i for Be targets, 
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Target 


Observable 


Geant4 


Penelope 


X 2 Ratio 






Factor 


Factor 


(G4/Pen) 


Be 


WT^(124keV) 


0.74 


1.10 


3.7 




$(124 keV) 


0.80 


1.09 


2.8 




r?(Si det.) 


0.84 


1.02 


3.1 




^(current int.) 


0.75 


0.91 


2.2 


Si 


1^^(124 keV) 


1.00 


1.08 


1.7 




|l(124keV) 


0.98 


1.10 


1.4 




r?(Si det.) 


0.98 


1.08 


1.3 




^(current int.) 


0.83 


0.91 


1.3 


Scint. 


^^(124 keV) 


1.07 


1.16 


1.8 




&(124 keV) 


0.94 


1.23 


4.8 




77(Si det.) 


0.92 


1.12 


0.74 




^(current int.) 


0.80 


0.97 


0.38 



ories which are improved versions of Moliere theory, and 
are therefore condensed algorithms requiring e.g. step 
sizes to be chosen very carefully |l2j |. Such algorithms 
generally perform adequately for small angle scattering. 
More recently, newer multiple scattering algorithms have 
become available, known as "mixed" algorithms which 
simulate hard collisions one by one (such as large-angle 
Mott scattering) and use a condensed algorithm to cal- 
culate the effects of soft collisions (such as small-angle 
Mott scattering and electron-electron scattering). Pene- 
lope belongs to the mixed class of simulation codes. 

We note that aspects of Penelope are included in the 
most recent versions of Geant4. However, the crucial as- 
pect of Penelope for the correct description of backscat- 
tering, which is the mixed algorithm multiple scattering 
code, is to date not included in Geant4. 



TABLE III: Overall scale factors and \ 2 ratios, comparing 
Geant4 to Penelope, under assumption of point-to-point un- 
certainty proportional to estimated normalization systematic 
uncertainty. 

where an absolute reduced \ 2 of 2.1 was seen for 188 
degrees of freedom. 

Overall normalization scale factors determined for 
both Penelope and Geant4 were generally found to agree 
within the normalization systematic uncertainties with 
unity, although deviations up to 25% are seen in certain 
cases. However, we note that, for observables which are 
not susceptible to extrapolation uncertainties, namely 
^-^^(124 keV) and ^(current int.), the Penelope scale 
factors are globally within 16% of unity. 

Penelope generally gives a lower x 2 than Geant4, in- 
dicating that the shape of the data is better described 
by Penelope. The exception is in the beam-energy de- 
pendence of T] for the scintillator data, for which Geant4 
gives a somewhat better fit. However, we believe this to 
be a coincidental cancellation, given that Geant4 gives a 
worse description of the data for the other two observ- 
ables for scintillator targets. 

The reason for the superior description of the data by 
Penelope is likely due to its treatment of multiple scatter- 
ing. In the case of backscattering, the multiple scatter- 
ing effects are dominated by large-angle scattering. The 
cross-section for large-angle scattering is dominated by 
Mott scattering. Multiple scattering algorithms, how- 
ever, do not necessarily include all collisions of a given 
particle. Such algorithms are referred to as "condensed" 
algorithms, where algorithms which include all collisions 
are referred to as "detailed" . Most particle physics simu- 
lation codes, such as Geant4, use multiple scattering the- 



VI. CONCLUSION 

Our new data on scintillator answer important ques- 
tions regarding the systematic uncertainties due to 
backscattering for a broad range of low-energy beta spec- 
troscopy experiments. Most importantly, the data agree 
well with models implemented in the Geant4 and Pene- 
lope Monte Carlo codes. 

Overall normalization scale factors were determined 
using a x 2 minimiziation technique. The resultant scale 
factors are generally found to agree within the normaliza- 
tion systematic uncertainties with unity. In some cases, 
discrepancies of up to 25% are seen. However, for ob- 
servables which are not susceptible to extrapolation un- 
certainties, the Penelope scale factors are always within 
16% of unity. In general, Penelope also gives lower x 2 val- 
ues than Geant4, indicating a better fit to the shape of 
the data in terms of energy and angle of the backscattered 
electrons, and in terms of beam energy. The reason for 
the superior description of backscattering by Penelope is 
likely due to its more accurate treatment of multiple scat- 
tering, employing a mixed algorithm treating large-angle 
scattering exactly, while using a condensed algorithm for 
small- angle scattering. 
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